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For the unitary operator, solution of the Schrodinger equation corresponding to a time-varying 
Hamiltonian, the relation between the Magnus and the product of exponentials expansions can be 
expressed in terms of a system of first order differential equations in the parameters of the two 
expansions. A method is proposed to compute such differential equations explicitly and in a closed 

■ form. 

o 
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^ ; I. INTRODUCTION 

in 

, For energy preserving finite dimensional closed quantum systems, the time evolution is well-known to 

^ ' be completely trivial. In fact, the linear and hermitian hamiltonian H admits the spectral decomposition 

. H = X^iLi with ei the energy eigenvalues corresponding to the stationary eigenstates and 

CO ' the solution of the associated Schrodinger equation in the \4>i{t)) is simply obtained by adding a phase 

, factor proportional to the corresponding energy eigenvalue: \(l>i{t)) = exp {—ieit/H) Even without 
resorting to the spectral resolution and to the change of basis that diagonalizes the Hamiltonian, the time 

, evolution of the quantum state \ip) can always be written as \ipit)) = U{t)\'ip{0)), with the unitary propagator 

■ U{t) — exp {—iHt/h) e SU{n) explicitly known. If Ai, . . . ,v4„, n = N"^ — 1, is a basis of su(A^), then it is 
\ possible to write U{t) in terms of some Eulcr-likc parameterization of SU{N). These parameterizations are 
. essentially ordered products of exponentials on SU(N) and for a Lie group there are as many such products 

I ' as there are Lie group decompositions. See |12, Q for a few examples of explicit choices on SU{N), 



X 



. |, I [TT[ |T| for SU{3) and g ||] for the related issue of how to uniformly parameterize the rotation groups. 

' Such decompositions are very useful for example in the generation of elementary gates in quantum computing 

^ I or more generally in the control of driven dynamics , but also in quantum state disentanglement [ p^ , 

<sJ '■ in the study of coherent states [l^ , in the solution of the Liouville-von Neumann equation and in the 

design of schemes for the numerical integration of differential equations on Lie groups |^ . 

The picture becomes more complicated when the Hamiltonian H is time-varying, because the energy (and 
thus the eigenstructure) is not conserved by the dynamics. The two expressions for the U (t) mentioned above 
also became time-dependent and go under the names of Magnus expansion [ p^ and Wei-Norman expansion 
p2| . The literature on the subject is quite vast, see and references therein, and deals especially with the 
convergence problem of the two series, which is not an issue here. 

We are interested here in studying how the two expansions relate, in particular how the parameters of the 
one series can be expressed as functions of the parameters of the other one. In nuce, such a transformation 
is already in the original papers of Wci-Norman, but its importance is made clear in p3| . In practice, it 
consists in studying a system of nonlinear differential equations having as variables the two sets of parameters 
contained in the two expansions. Such a system of differential equations is called the Wei- Norman formula 
and it corresponds to the Jacobian of the change of coordinates i.e. of the transformation from single 
exponential to product of exponentials. As such, it is linear in one of the two sets of coordinates although 
it depends on the point of application, expressed nonlinearly in terms of the other set of parameters. 

The Wei-Norman formula appears in several different contexts in the literature, see |[ |^, |l^ just to 
mention a few. In this study, we propose a method to compute it explicitly and systematically for any 
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dimension, based only on the structure constants of the Lie algebra. Such method seems to be new. The 
rationale behind it is a technique to compute in closed form one parameter groups of automorphisms i.e. 
exponentials of the matrices of the adjoint representation of any linear Lie algebra. Such technique is adapted 
from Ip^ . As an example, we compute two Wei-Norman formulae for su(2) for the same basis obtained from 
the Pauli matrices but for different ordering of the basis elements (i.e. different choices of Euler angles) and 
a formula for su(3). As we will see, the explicit symbolic computations rapidly grow with the dimension N; 
however numerical simulations of the dynamics can always be carried out easily. 



II. MAGNUS EXPANSION VERSUS PRODUCT OF EXPONENTIALS EXPANSION 



Consider a quantum mechanical system with Hamiltonian H{t) which is a continuous function of time. 
Assume H{t) is Hermitian and belongs for all t to the finite dimensional Lie algebra 5u(A''). If we choose 
a basis of skew-symmetric matrices Ai, . . . An for su(A^), then [Ai, Aj] = c'^jA^^ where c^j are the structure 
constants of su(A^). In the following we will assume that the summation convention over repeated indexes 
is imposed only for the greek letters, not for the latin. Since Ai, . . . ,An are time independent operators in 
su(A^), H{t) can be written in this basis as H{t) = u^{t)Afj_ with analytic functions of time. Let U{t) 
be the time evolution operator of the Schrodinger equation: 

^l\^it)) ^ H{t)m)) \m)^m 

Then = U{t)\'ipo) where U{t) in the Magnus expansion is given by the formal expression 

U{t) = Texp u''{T)A^dT 



(1) 



where T is the Dyson operator and exp, the exponential map for SU (N), is the ordinary matrix exponential. 
The Wei-Norman formula relates (|l|) with the expansion as a product of exponentials, i.e. it affirms that 
(0) can be written locally around the identity of SU {N) as 



Uit) - exp {^\t)A,) ... exp (7"(i)^n) 
In other words, it is used to relate the flow of the solution of the differential equation 

. d 



(2) 



i-U{t) = H{t)U{t) 



U{0) - / 



with the product of one-parameter subgroups 7* lifted to SU {N) and used as a coordinate chart on the 
group itself. If, in the original formulation of Wei-Norman, the same set of infinitesimal generators as the 
basis is chosen, nobody forbids to alter the order of the Ai in or even to consider any alternative 
fundamental representation of the group, for example one of the Euler-like parameterizations introduced 
in m |ri|, ^Sj. We will stick to (||) for notational convenience, but develop the corresponding of the 
Wei-Norman formula for the Euler ZYZ parameterization for the SU (2) case in the example of Section ^ 

The Wei- Norman formula consists in expressing the functions 7'(i) in terms of the u^{t) via a system of 
differential equations 



S(7\...,7") 



7X0) = 



(3) 



with the n X n matrix ^ analytic in the variables 7*. Since the adjoint maps can all be written as 



(4) 



1=0 
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the matrix S of elements (S)fei = is defined in terms of the 7' and of the structure constants as: 

m 



(5) 



When the Ai form a basis of su(A^), the matrix S assumes also the meaning of map between canonical 
coordinates of the first kind and canonical coordinates of the second kind (^), see ||2l[]. In this case, since 
7*(0) = 0, S(0) = / and thus S is locally invertible. However, because of the semisimplicity of SU{N), all 
parameterizations lead to a Wei-Norman formula that is subject to singularities and as such has only 
a local validity. By inverting S, equation (j3|) assumes the more traditional aspect of a system of first order 
differential equations in the 7' variables: 



7" 



s(7\...,7"r' 



7'(0) = 



(6) 



If the time evolution of one of the two vectors of coordinates 7* or it* is known, the formula3(||) or (||) can 
be used to obtain the other one. While (||) is global, (|^) is valid only as long as det(S) 7^ and thus the 
nonsingularity of S needs to be checked at the point of application. Another weak point of the Wei-Norman 
formula is that the system of differential equations is nonlinear in the 

Finally, (|4|) as it stands involves infinite sums. In the following paragraph we provide a method for 
transforming these infinite series into finite sums. 



III. COMPUTATION OF THE EXPONENTIALS IN THE ADJOINT REPRESENTATION 

The adjoint representation of a linear Lie algebra like su(iV) is a derivation of the algebra and as such it 
is the infinitesimal generator of a one-parameter group of automorphisms. In the basis Ai, . . . , An of the Lie 
algebra, the one-parameter groups of automorphisms are described by the matrix exponentials e'^'^^i |^ . 

There exist several techniques to compute the exponential of a square matrix. The simplest is perhaps the 
use of the Sylvester formula [Q, p. 233, which requires only the knowledge of the eigenvalues of the matrix. 
Alternatively, one can use the spectral decomposition of normal operators i.e. compute the Jordan form by 
similarity transformations, but this requires the knowledge of eigenvalues and eigenvectors. The method we 
use here is taken from [Q. It relies on the Cayley-Hamilton theorem and consists in expressing the series 
expansion of e'^'^^i in terms of the first n — 1 powers of ad^^ with suitable coefficients depending on the 
coefficients of the characteristic polynomial of ad^. and on 7^. Of the different methods that can be used 
to express a matrix exponential, in our opinion this is the most suited for the adjoint representation, as the 
powers of ad^i^ are immediately expressed in terms of the structure constants of the Lie algebra. In fact, if 
c'^j are the structure constants of the Lie algebra, the matrices corresponding to the basis elements Ai are 
adyi; — Mi of elements {Mi)kj = c^^ and the l-th power of ad-Ai is given by 

Af/ = (Ml) , . = cf „ cff], . . . c't'-' cfr' (7) 



Any matrix B G 5u{N) can be written as i? = If we identify B with the n-dimensional coordinate 

vector B — ... 6"]^, then [A^, B] ~ ad^^^ — {Mi)b, i.e. the Lie bracket gives another column 
vector 

while the Lie bracket oi D ~ d = [d} ... d"] ^ with B looks like: 

[D, B] ^ adob - d^iM.jb = [d'^cl^b^^ d'^cl^b'^ . . . d•^c';^,fi^^]'^ 
As an example, compute [B, Ai] — ~[Ai, B]. Since Ai ~ e^, the standard basis vector of M", we have: 

[B, A,] ~ b''{M,)k^ = b"- [4,; . . . c:,] [cl . . . c^] ^b'' = -{M,)k,b'' ~ B] (8) 
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A. Exponential of the adA matrix 

Since the adjoint representation corresponds to a derivation of the Lie algebra, given A e 5u{N) 

7 e M 



ad^ = -7- (Ade^A) 



7=0 

and ad^ = M is the infinitesimal generator of the corresponding one-parameter group of automorphisms 

00 00 ^ 

Ad,A. = e^^^- = e^^ = 5^ ^ad^ = J2 (9) 

fc=0 ' fc=0 

To obtain an explicit closed form expression for , we use the method of which is based on the Cayley- 
Hamilton theorem i.e. on expressing the exponential of a n x n matrix in terms of its first n — 1 powers. In 
fact, if M satisfies its own characteristic equation 

M" = aal + aiM + aaAf^ + . . . + a„_iAf"-i (10) 

then (m can always be written as 



ef^'i ^Y.PkM" (11) 



for suitable Pk — Pkio,Oi ■ ■ ■ ,o-n~-i, 7)- In p^ , the values Pk are obtained via a Laplace transform method 
that is briefly recapitulated here. From (O) and (foh, one can write /3„_i as the infinite series 



Pu-i = E 4"^- (12) 
j=o ^■ 

with aj, j — n ~ 1, n, n + 1, . . . satisfying the linear recurrence relation with constant coefficients 

n-l 

Oij+i = ^ a„_fe_iaj_fc = aoaj-„+i + aiaj_„+2 + • • • + fln-iaj (13) 

(14) 



*;=0 

Jo ifj<n-l 
1 1 a j = n-l 



If we Laplace transform termwise /3„-i, £ |^q;j| — ■jjtt, then the infinite sum ( p^ becomes a 2-transform 



m the s variable, Z{s) = '^JLo '^j^ ^-^^ 



Multiplying both sides of (|l3| ) by s summing for j from n — 1 to 00 and completing the series to zero 

by means of the initial conditions (114) , we obtain 



j=0 j=0 



from which, if s^, i = 1, . . .r < n, are the eigenvalues of det(s/ — M) of multiplicity m^, the difference 
equation (p^-(p^ can be written in terms of this transform as 
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or, expanding in partial fractions 



where 



r rui 



^ is- SiY 



[nil - j)\ \ds 



(z(s)(.-..r) 



Thus antitransforming 



(15) 



4=1 



with Cii^j) = X^j^i C'ij ■ The Pki k = 0, 1, . . . , n — 2, are then obtained by differentiation of /3„_i with 
respect to 7: 



-fc-i 



-k-l 



(16) 



B. Exponentials of the basis elements adA^ 

For the basis elements Ai € 5\x{N), the powers of the matrices Mi = adyi- are computed in terms of 
the structure constants in (0). Thus, using the notation /3q /jj^'L 1 for the coefficients obtained in 
(p^)-(|l6|) relative to the matrix M^, we have 

n-1 

^ Mr '^i Ml S Ai2 • ■ ■ J 
r=0 

where it is intended that cj^,, cf,^, . . . cfp^ = 5'^ for r = and „ d!'], . . . cfj"^ — c% for r = 1 (the lower 

index in gives the number of times the structure constants c*^ appear in the corresponding term). Also 

the product of matrices '^'^^^ and can be expressed in terms of the /3q ' , . . . , P^n-n /^o ' j ■ • ■ i P^n-i and 

of the structure constants as follows: 

^.^ad^^^Vad,, VM.^.'mA ^ ^ /3M/3i'lctiCr;,,...cfr^cr,,cr^,...cr;- (18) 

n 

r, s=U 

with a similar convention as above for r, s = 0, 1. Notice how also the "concatenation" of c*, and C;*^ respects 
the summation convention (indeed ( |l8| ) represents an ordinary product of square matrices). 



IV. EXPLICIT EXPRESSION OF THE WEI-NORMAN FORMULA 



With the notation introduced in Section [I], the matrix ^ of (H) is given by 



ei e^'^d^iea ... UZi e-^^'^^^^^r. 

With the closed form expression (|l^) for ^'^-^i , it is possible to compute S explicitly without any infinite 
summation. Looking at (||) and (|l^) it is clear what we have to do: from the matrix products of the '^^^i 
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compute all the column vectors corresponding to Y[i=i e^^'^'^'^^Aj, j = 2, . . . ,n, in the same fashion as (^) 
and then regroup them along the n basis directions Ak . This can be done explicitly. 



. . c 



r=0 
n-1 



■-2 3 



I"! ,'"2=0 



(19) 



n-l 



n-1 



n 7'adA, 4 _ nil] fjln-l] k Mri-l ^ /'""-i"! 

i— 1 ri,...,rTi-i=0 

The only free index in (|l9|) is fc i.e. each Hti '^'^^^ n-dimensional vector as expected. The 

expression of the n x n matrix S one obtains by stacking together the columns vectors of ([l9|) is 



■-, spn-l M l Mr-l 
^ Z^r=0 '-iMl ■ ■ ■ ^Ij 



Z^n r„„i=0 ^'"1 ■ ■ ■ Cj^^ ... 



Z^ri,...,r„_i=0 '-'''1 ■ ■ ■ Pr„-i ■ ■ ■ <^n-l 



(20) 



V. EXAMPLE: su(2) 
A skew-symmetric basis for su(2) is obtained from the Pauli matrices 





1 









'1 " 


(Tl = 


1 




i 




-1 



for example by taking Aj — ^cjj, j — 1, 2, 3, i.e 



Ai = 



1 



i 

1 



A. - 



1 

-1 



As = 



i 

-i 



(21) 



and it corresponds to all real structure constants — C23 — c^i — 1- The corresponding adjoint matrices 
are 





"0 


" 




" 


1" 




"0 


-1 0" 


ad^i = 





-1 


ad^2 








adA3 = 


1 










1 




-1 














whose exponentials are already known from the literature to be 

7^adAi 



1 

cos(7^) — sin(7^) 
sin(7^) cos(7^) 



o7 ad^ 



cos(7^) sin(7^) 

1 
-sin(7^) cos(7^) 



cos(7'^) — sin(7'^) 
sin(7'^) cos(7'^) 
1 



The Magnus expansion (|l]) is given by U{t) = T e 



dr (with = u*(t)). 



A. Wei-Norman formula for canonical coordinates of the second kind 

In the product of exponentials, choosing the order given by the cardinality of the index gives the canonical 
coordinates of the second kind on SU{2): U{t) = '^^e* "^^e^ (again with 7* = 7*(0)- For this choice, 
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the Wei-Norman formula 



reads as: 



1 sin(72) 

cos(7^) — cos(7^) sin(7^) 

sin(7-^) cos(7^) cos(7^) 



whose inverse can be computed explicitly: 



1 sin(7-'-) tan(7^) — cos(7-'-) tan(7^) 
cos(7^) sin(7^) 
— sec(7^) sin(7^) cos(7^) sec(7^) 



From 



the determinant of ^ is simply 



det ^ = cos 7 



(22) 



(23) 



and thus the singularities of the representation are 7^ — tt/2 + kn, fc e Z. While (H) (and (p2[)) is valid 
everywhere, in the singular points the formula (^ cannot be applied (i.e. ( p3| ) is not defined). 



B. Wei-Norman formula for the ZYZ Euler angles 



When expressed in the basis (|2l|), the ZYZ Euler angles correspond to the product of exponentials U{t) ~ 



this case to 



e ' '"■'e^ (compare with the expression of Section V A ) . The Wei-Norman formula corresponds in 



ea e'^''^<^'^3e2 



63 



— sin(7-^) cos(7^) sin(7^) 

cos(7^) sin(7^) sin(7^) 

1 cos(72) 



and its inverse 



Since 



— cos(7-^) cot(7^) — sin(7-^) cot(7^) 

— sin(7-^) cos(7^) 
cos(7^) csc(7^) sin(7-^) csc(7^) 



det(^) = sin 7^ 



the singularity has now moved to 7^ — kn, fc g Z, as is well-known for such a parameterization. Thus 
can be used everywhere except in the identity U (0) = /. It is worth emphasizing that it is a fundamental 
topological fact that singularities cannot be avoided in a minimal parameterization of a semisimple Lie 
group. One possible way to get around the problem is obviously to use "redundant" parameterizations like 
quaternions. 



VI. EXAMPLE: su(3) 

The following skew-Hermitian basis obtained for example by considering su(3) as the compact real form 
of the classical A2 algebra with its Cartan basis is more suited for the analysis of three level systems than 
the Gell-Mann basis normally used in particle physics. 





I 














"000" 




" 1 


0" 




"0 


i 





iHi = 





— i 







IH2 = 


i 


X12 — 


-1 





Y12 = 


i 

























































1" 






"0 i 




"0 


0" 




"0 





0' 

















Yi3 = 







X2i = 





1 


Y23 = 








i 








-1 












i 




-1 










i 






(24) 
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The nonnuU structure constants c?,- — — c*^, are: 

*j J* 

'-IS ~ ^-41 ~ '-34 ~ ^' '-IS ^ '-ei ~ '-56 ^ ^' '-27 ~ ''82 ~ '-78 ~ ^ 
C45 = ^84 — ^58 — 1; '^46 = '^74 — ^67 — "1; ^35 " '^73 — 4,7 " '^36 = ^83 " '^68 — ~^ (^5) 

'^17 — '^81 — "1; '^78 ~ '^23 = '^42 = ^34 = 0; ^25 = '^62 ~ 1; ^56 " ^ 

Note that because of the last row of (^ ) the presentation of ad^ is not skew-symmetric, hke it would have 
been had we chosen the Gell-Mann basis. 

For su(3), the basis of the adjoint representation is 





"0 




















" 




"0 




















" 

































































„3 

























'-24 














Ml = 










'-13 














„5 
'^16 




















'-23 














„5 
'-26 
























^15 

























'-25 



























































'-28 






















„8 
Cl7 


. 






















„8 
^27 


. 



" 








'-34 























































'-32 



































C37 



































C35 























„8 
'-36 












'-43 














" 


























'-41 


'-42 

































































C48 




















C47 




















C46 




















„8 

C45 














" 














''56 





" 




" 











'-65 








" 



















'-56 






















'-65 































„3 

'-57 




























'-68 


M5 = 






























'-58 








„5 
^61 




„5 
'-62 


















'-67 










'-51 


'-52 






















































'-53 




























'-64 
















. 








^54 






















'-63 


















' 




















" 























2 
'^78 














„3 

C75 


























'-76 














C73 


























„6 

C74 






































„8 

LC71 


„8 

C72 





















" 




















0" 




















^2 
^^87 




















„3 
'-86 




















'-85 




















„5 
'-84 




















'-83 

















'-81 


'-82 













































The characteristic polynomials and corresponding eigenvalues of the basis elements Mi 



Aet{sl- Mi) = 




95*^ + 24s'* + 16s2 
Ss*^ + 13s4 + 6s2 



for i = 1 
for i = 2, 3, 4, 7, 
for i = 5, 6 



(26) 
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{0, 0, ±i, ±2i, ±2i} for i = 1 
eig(Af,) = { {0, 0, ±i, ±i, ±2i} for i = 2, 3, 4, 7, 8 
{0, 0, ±i, ±i, ±y/6i} for i = 5, 6 



Inserting the coefficients of the characteristic polynomials (| 
following values for the , fc = 0, . . . , 7 and i = 1, . . . , 8: 

1 



into (p^, one gets from (|5|) and ( |l^ ) the 



[i] 



4" 



4'' 



for i = 1 , . . . , 8 

7* fori = l,...,8 

' ^ (54 - 64 cos(7') + 10 cos(2 + 3 7* sin(2 7*)) 
< 5^ (81 - 80 cos(7'') - cos2(7*) - 247' sin(7') + sin2(7')) 
jig (325 - 324 cos(70 - cos(y6 70 - 90 7' sin(7*)) 

^ (2167^ - 67* cos(2 70 - 256 sin(7*) + 23 sin(2 7*)) 
^ (81 7' + 247' cos(7*) - 104 sin(7') - cos(7*) sin(7^)) 

137' I 37* cos{7^) 69 sin(7^) sln(^/6^^) 
~6~ ^ 5 25 150 

(81 + 47 cos(2 70 +cos(70 (-128 + 307'' sin(7'))) 
26 cos(7'') - cos2(7') - 15 7' sin(7') + sin2(7')) 
(200 - 198 cos(7'') - 2 cos(V670 - 105 7' sin(7')) 

gig (3247' - 30 7* cos(2 7') - 512 sin(7') + 109 sin(2 7')) 
^ (277* + 157' cos(7*) - 41 sin(7*) - cos(70 sin(7*)) 
(600 7' + 3157' cos(7') - 909 sin(7') - ^/6 sin(y6 7')) 

^(9 + 7cos(2 7*)+cos(7*) (-16 + 67* sin(7'))) 
^ (9 - 8 cos(7') - cos2(7') - 67' sin(7') + sin2(7*)) 
jig (25 - 24 cos(7') - cos(\/6 7*) - 15 7' sin(70) 

' gig (367' - 67* cos(2 7*) - 64 sin(7') + 17 sin(2 7')) 
^ (97* + 67' cos(7') — 14 sin(7') — cos(7') sin(7')) 



144 

^(27 



gig (150 7' + 90 7' cos(7') - 234 sin(7') - \/6 sin(V6 7')) 
and, from (p+|), the exponentials are: 



p7iadAi _ 



10 
1 

cos(27i) 

sin(27i) 
















-2cos(7^) sin(7^) 
cos(27i) 













cos(27i) 
sin(27i) 







for 1 = 1 

for i = 2, 3, 4, 7, 8 
for i = 5, 6 

for i = 1 

for i = 2, 3, 4, 7, 8 
for i = 5, 6 

for 1 = 1 

for i = 2, 3, 4, 7, 8 
for i = 5, 6 

for i = 1 

for i = 2, 3, 4, 7, 8 
for i = 5, 6 

for i = 1 

for i = 2, 3, 4, 7, 8 
for I = 5, 6 



for i 
for i 
for i 



1 

2, 3, 4, 7, 8 
5, 6 








~2 cos(7-'-) sin(7-'-) 
cos(27i) 













cos(7^] 










sin(7^) 



■sin(7-^) cos(7^). 



"1 

10 

cos(72) sin(72) 

-sin(72) cos(72) 

cos(72) 

sin(72) 















-sin(72) 

cos(72) 



cos(272) -sin(272) 
sin(272) cos(272) 
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o7 adAo 



cos(273) sin^(73) sin(273) 

10 

10 

-sin(27^) cos(7^)sin(7^) cos(27^) 









COS (7^) 


-sin(7^) 











cos(7^) 


- sin (7^) 








sin(7^) 


cos (7^) 




■ 






sin (7^) 


cos(7'^). 



■cos(274) 


sin(274) 








sin2(74) 
1 

■ cos(7^) sin(7'') 








- sin(274) 


cos(274) 
















cos(7"^) 



sin(74) 










COS (7^) 










sin(7^) 



- sin(7'') cos(7^) 










— sin(7**) 



cos(7^) 



1+2 cos(^/67^'^) -1+005(^6- 
3 3 
2 (-l+cos(V67^)) 
3 











2+cos(\/67°) 
3 






sin(\/6 7" 






cos (7^) 





sin(7^) 









cos(7^) 






- sin(7^) 



7" J 

sin(V675) 





cos(-\/67^) 









— sin(7^) 





COS (7^) 









sin(7^) 






COS (7^) 



1+2 cos(\/67'') 
3 

2 (-l+cos(V6 
3 



^ sin(V67^) 






-l+cos(\/67'') 
3 

2+008(^67") 
3 





sin( V6 7*^) 
V6 










COS (7®) 








sin (7^) 









cos(7^) 



sin (7^) 




sin(V67'')) 

sin(-\/67®) 



cos(-\/67^) 













— sin(7^) 



cos(7^) 







— sin(7^) 






COS (7^) 



g7''adA7 



1 

sin2(77) 










cos(27'^) 








.cos(7'^) sin(7'^) — sin(27'^) 






cos(7^) 


— sin (7'^) 











cos(7'') 


— sin(7'^) 








sin(7'^) 


003(7'^) 











sin(7'') 



COS (7'^) 

1 




sin(27'^) 



cos(27'^). 
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o7 adAg 



1 

sin2(78) 









003(278) 







-003(78)3111(78) 3111(278) 







003(78) 



— 3111(78) 









003(78) 






3111(78) 

- 3111(78) 



0" 
3111(278) 



3111(78) 003(78) 



003(78) 









003(2 78) 



Finally, the expre3sion for the Wel-Norman formula Is: 



10 
1 

003(27! 



3ln(27i - 7^) 











sln(273) 


- 3111(27! - 7^) 003(27^) 
003(27! - 7^) 003(27^) 






















$17 


$18- 





y|sln(^/675) 





$28 





$36 


$37 


$38 





$46 


$47 


$48 


^55 


$56 


$57 


$58 


^65 


$66 


$67 


$68 


^75 


$76 


$77 


$78 


^85 


$86 


$87 


$88. 



where the rather oumbersome expliolt expression of various terms Is Included below for sake of oompleteness. 

$55 = 003(27! + 7^) 003(7^) 003(7^) - sln(27! + 7^) 3ln(7^) sln(7'') 

^g5 = 003(7^) 003(7^) sln(27! + 7^) + 003(27! + 7^) 3ln(7^) sln(7**) 

$75 = — 003(7! — 27^) 003(7^) sln(7'^) + 003(7"^) 3ln(7! — 27^) sm(7**) 

$85 = 003(7"!) 3ln(7! — 27^) sm(7^) + 003(7! — 27^) 003(7*^) sm(7'!) 



$16 = 

$36 = 

$46 = 

$56 = 

$56 = 

$76 = 

$86 = 



sm 



(V67'' 



2V6 

sln(\/67^) 



4^/6 



(2 + 003(2(73 - 7^)) + 003(2(73 + 7^))) 

(- 003(27! - 7^ + 273 - 27^) + 003(27! - 7^ - 27^ + 27*) + 003(27! - 7^ - 2(7^ + 7*))- 



- 003(27! - 7^ + 2(73 + 7*)) +4003(27! - 7^) sln(27'')) 
sin(-v/67^) 



2\/6 



(003(27! - 72 - 27"!) - 003(27! - 7^ + 27"^) - 2003(27! - 7^) 003(27^) sln(273)) 



- 003(^67^) (003(7^) 003(7"^) 3ln(27! + 7^) + 003(27! + 7^) 3ln(73) 3ln(7'^)) 
0O3(\/67^) (003(27! + 7^) 003(7^) 003(7"^) - sln(27! + 7^) 3ln(73) sln(7"!)) 

- oo3(\/67^)(oos(7'*) sln(7! - 27^) 3ln(73) + 003(7! - 27^) 003(7^) 3ln(7^)) 
oos(-\/67^)(— 003(7! — 27^) 003(7"*) sln(73) + 003(7^) sln(7! — 27^) 3ln(7'')) 
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^17 = cos(27^) 003(76) sin(27*) 3111(7^) - 003(7^) sin(273) 3111(76) 
^37 = - 003(27^-72) 003(27^)003(76) 3111(7^)+ 3111(27! -7^)- 

• (003(76) 3111(27^) slii(27*) sln(75) + 003(27^) 003(7^) sln(76)) 
^47 = - 003(27*) 003(76) 310(27! - 7^) 3111(7^) - 003(27! - 7^) • 

• (003(76) 3111(27^) sln(27*) 3111(7^) + 003(27^) 003(7^) 3111(76)) 

^57 = - 3111(27! + 7^) (003(7^) 003(7^) 003(76) sln(7*) + 003(7^) sln(7^) 3ln(7^) 3111(76)) + 

+ 003(27! + 7^) (003(7*) 003(7^) 003(76) sln(7^) — 003(7^) 3ln(7*) 3111(7^) 3111(76)) 
^67 = C03(7*) sln(7^) (003(7^) 003(76) 3111(27! + 7^) + 003(27! + 7^) 310(7'"^) 310(76)) + 

+ 003(7-'') sin(7*) (003(27! + 7^) 003(7'"^) 003(76) - 3111(27! + 7^) 3111(7^) 3111(76)) 
= 8111(7! - 27^) (- 003(7^) 003(76) sln(7-') 3ln(7*) + 003(7^) 003(7*) 3111(7'"') 810(76)) + 

+ 003(7! - 27^) (003(7^) 003(7*) 003(7^) 003(76) + 3m(7^) 3ln(7*) 310(7^) 3111(76)) 
^87 = 003(7^) 003(7*) (- 003(7^) 003(76) 310(7! - 27^) + 003(7! - 27^) 3111(7^) 3111(76)) - 

- 3ln(7^) 3ln(7*) (003(7! - 27^) 003(7^) 003(76) + 3111(7! - 27^) 310(7^) 3111(76)) 

^81 = C03(27^) (003(76) sln(27^) 3ln(7^) + 003(27^) 003(7^) sin(27*) 3111(76)) + 

+ ^ ((2 + 0O3(v^75) 003(^676)) sln2(73) + 

+ 003(27^) (003^(7*) (-1 + 003(^67^)003(7676)) +3 3ln2(7*))) 3ln(27^) 
= 5^!^(2 + oos(x/67^) oo3(^^7')) 

^83 = - 003(27'^) (003(27! - 7^) 003(27*) 003(7^) 3111(76) + 

+ 3111(27! - 7^) (003(27^) 003(76) 3111(7^) - 003(7^) 3ln(27^) 3ln(27*) 3111(76)) ) 
^sin(|/) ^ co3(^^75) cos{Vh')) ■ 

■ (003(2 7*) 3ln(2 7! - 7^) 3ln(2 7^) + oo3(2 7! - 7^) 3ln(2 7*)) 

^84 = 003(27'^) (- 003(27*) 003(7^) 3111(27! - 7^) 3111(76) + 

+ 003(27! - 7^) (003(27^) 003(76) 3ln(7^) - 003(7^) sln(27^) 3ln(27*) 3111(76))) 

_^sin(27^ (-4 + 003(767^) 003(7676)) (- 003(27! - 7^) 003(27*) 310(27^) + 3111(27! - 7^) 3ln(27*)) 

^85 = 003(27"^) (- cos(7*) 3ln(73) (oo3(7-'^) 003(76) 3ln(27! + 7^) + 003(27! + 7^) 3ln(7-"^) 3111(76)) + 
+ 003(7^) 3ln(7*) (- 003(27! + 7^) oo3(7''^) 003(76) + sln(27! + 7^) slii(7''^) 3in(76))) + 

_^sm{^_)_ (003(7676) (003(7^) 003(7*) §111(27! + 7^) + 003(27! + 7^) 3111(7^) 3ln(7*)) 3111(767^) + 
76 ^ 

+ (003(27! + 7^) 003(7^) 003(7*) - 3111(27! + 7^) sln(7^) 3ln(7*)) 3111(7676)) 

^86 = C03(27^) (- 3ln(27! + 7^) (003(7^) 003(7"') (-03(76) 3ln(7*) + 003(7*) 3111(7^) 3ln(7-'^) 3111(76)) + 
+ 003(27! + 7^) (003(7*) 003(7^) 003(76) sln(7^) — 003(7^) 3ln(7*) sln(75) 3111(76))) + 

_^sm(^j_)_ (003(7676) (- 003(27! + 7^) 003(7^) oos(7*) + 3111(27! + 7^) sln(7^) 3ln(7*)) sln(767^)+ 

76 ^ 

+ (003(7^) 003(7*) 3111(27! + 7^) + 003(27! + 7^) 3ln(7^) 3ln(7*)) 310(7676) J 
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^88 = cos(27^) (sin(7^ - 27^) (- 003(7^) cos(7*^) sin(7^) sin(7^) + 003(7^) 003(7^) sin(7^) sin(7^)) + 
+ cos(7^ - 27^) (cos(7^) 003(7"*) 003(7^) 003(7'') + sin(7^) 3in(7'*) sin(7^) sin(7''))) 



This time the singularities of the Wei-Norman formula are obtained from the zeros of the following function: 



To be able to describe a time varying dynamics in terms of simple unitary operations is an important issue 
in quantum mechanics and it is foreseen that it will be a crucial one in quantum computation. The method 
we propose here relates the one parameter flow of the Schrodinger equation with an arbitrary decomposition 
of SU (N) by computing the Jacobian of the coordinate transformation. It is worth emphasizing that if 
the symbolic expression of the Wei-Norman formula rapidly explodes with the dimension of the system, its 
numerical integration can be easily and efficiently handled. Furthermore, because of the way it is formulated 
(each parameter in (|l|) and (|^) lives on M* or on §*), the numerical integration will give a solution which 
respects the group structure, regardless even of the rounding error. Finally, notice that, as a byproduct, the 
explicit expressions for the exponentials of the adjoint operators should be useful in the study of the solution 
of the density operator equation, which is often expanded in terms of unitary superoperators. 
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